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1. Introduction 

Phenomenological studies of lipid bilayer membranes have become imperative, 
for a large amount of experiments reveal that the matrix of biomembrane plays an 
important role in some non-linear phenomena such as shape fluctuations of cells 
and oscillations of the membrane potential in an external electric field (1-3). In 
approaching the problem, however, one is immediately confronted with complexities 
because of the concurrence of the elastic- viscous motion and the change in shape of 
the material. The lack of a theoretical framework tailored for non-Newtonian fluid 
membranes impedes progress in this research area, though the general geometrical 
method in continuum mechanics has been well established. The purpose of this 
t^j- ! paper is to present such a framework. 

The membranes constituted by two layers of rod-shaped molecules are identified 
with hydrotropic liquid crystals in the laminar phase. From this point of view, one 
assumes in the present paper that the dynamic theory of liquid crystals developed 
during the last 30 years is applicable to the membranes. In fact, the experimental 
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CN . observations mentioned above would not be surprising if one thinks of the anisotropy 
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feature of the material, which, as is well known, leads to instabilities in the liquid 



f — ; crystal bulk. 

Due to the pioneering research of Oseen [4], the fundamental contribution of 
Frank [5] and the important supplementary work of Nehring and Saupe [6] , a cur- 
vature elasticity theory of liquid crystals has been established. It is essentially 
a phenomenological theory. A dimensionless unit vector named "director" is de- 
fined to indicate the preferred orientation of the molecules. The energy stored in 
the material during elastic deformation is thus expressed in terms of this vector. 
Extending the theory to thin films in particular yields the well-known curvature 
elasticity equation of membranes [7] which has been used in fruitful studies on the 
hydrostatic form of lipid vesicles [8] . 

The dynamic theory of liquid crystals [9-13], mainly uniaxial crystals, stems 
from such a mechanical base. It assumes rotation of the molecules, indicated by the 
director, superposed on the mean flow of the molecules themselves. The theory was 
used to investigate instabilities of liquid crystals and has interpreted successfully 
the Williams' domain structure and the homoetropic textures as well in nematics 
[14-22]. 

One adapts the theory to the particular context by neglecting non-uniformities 
across the membrane thickness and replacing the three-dimensional film by a two- 
dimensional surface. It is exactly the method used in the theory of thin shells 
[23]. The simplification is acceptable in most cases of biocells, for the membrane 
thickness (40-80 A) is sufficiently small with regard to the cell diameter (1-10 /im). 
Moreover, the tilt angle of the long molecular axes with respect to the vector normal 
to the surface is supposed to be small and hence the elasticity theory is restricted 
within the extent of linear approximation. With these considerations, one obtains 
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a particular version of the dynamic equations, in which geometric invariants such 
as the fundamental magnitudes and the curvatures are used to describe the shape 
of the membrane, while the displacement of the configuration of the membrane is 
treated separately from the internal movement of the particles inside the membrane. 
The attempt to construct such a group of equations, so far as one knows, was made 
10 years ago by Peterson [24]. He established, however, only the dynamic equations 
of a solid film and the hydrodynamic equations of a Newtonian fluid membrane, 
separately. This left a gap between these two limiting cases. One will see that by 
introducing the theory of liquid crystals the gap is naturally bridged. 

Recent tether experiments support the opinion that slip between the two molec- 
ular layers causes inter layer viscous dissipation [25] which might not be negligible 
in dynamical studies of bilayer membranes [26]. One recognizes that this kind of 
internal dissipation is excluded from the present theoretical model. 

This paper is the first part of a series devoted to present the general electro- 
hydrodynamic equations of polar viscoelastic fluid membranes isolated from their 
surroundings. The hydrodynamic theory of liquid crystals will be recapitulated in 
the next section, then it will be adapted to the context of bilayer membranes in the 
third section. The connection of such a membrane to the liquid surroundings will 
be considered in subsequent papers. 

2. Outline of the Ericksen-Leslie theory 

The hydrodynamics of liquid crystals was primarily established by Ericksen 
[9]. He assumed an intrinsic motion of the non-spherical molecules and used a 
time-dependent director to indicate it. Because of the material symmetry, he con- 
structed the constitutive relationships of the dynamic variables involved and hence 
interpreted explicitly the dynamic conservation laws for uniaxial liquid crystals. 
Leslie [10] was able to take into account the antisymmetric stress of the anisotropic 
continuum [27,28], which had been overlooked in the expressions of Ericksen, and 
thus provided us a refined version of the theory. 
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2.1. Mass balance 

The principle of mass conservation holds in the convected volume R 



~dt 



JJJpdR = (1) 



R 

where p is the mass density and t is the time. 

2.2. Conservation of momentum 

Let v denote the barycenter velocity of the molecules. Conservation of linear 
momentum in the motion of the continuum is assumed, which asserts 

iJJJ? vdR = JJJ FdR+ JJ ds ' T (2) 

R R s 

where F is the body force, T the stress tensor and s is the surface enclosing the 
considered volume. 

The body force consists of the gravitational attraction and electromagnetic 
forces [9] 

F = F 5 + F e ' m (3) 
The gravitational force is naturally written as 

F 9 = pg (4) 

where g is the gravitational acceleration. 

The force generated in an external magnetic field is given as [9] 

F m = (M-V)H (5) 

where H is the external magnetic field, M is the magnetization and V is the vector 
derivative operator. 

The most general relationship of M linear with the magnetic field is [9] 

M = X ±H + (x|| - X ±) (d H) d = X ±H + X a (d • H) d (6) 

where d is the director, x\\ an d are the constant magnetic susceptiblities in the 
directions parallel and perpendicular to the director, respectively, and Xa is the 
anisotropy. 

The force generated by the electric field is not expressed explicitly in the 
Ericksen-Leslie theory. In electrohydrodynamic studies of liquid crystals, however, 
the Lorentz force is commonly considered as a cause of instabilities [17-22]. We do 
not go far towards calculating completely the force in the electric field, but, inaccu- 
rately, equate it to the Lorentz force, which amounts partly to the ponderomotive 
force in liquid dielectrics [29]. That is 

F e = p e E (7) 
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where p e is the density of injected electric charges and E is the electric field. 
The stress tensor consists of three parts [9] 



T = T l + T e + T v (8) 

where T l is the isotropic part of the stress, T e is the stress arising from elastic 
deformation and T v is that from viscous motion. 
The isotropic part is expressed as 

T l = -al (9) 

where a is a constant and I is the unit tensor. 

The elastic stress may be deduced from the density of reversible deformation 
work W d [9], i.e. 

dW d 

T ^-< Vd »'9(dV) < 10 > 
It is suggested in the Oseen-Frank-Nehring-Saupe theory that 

W d =h (V • d) - k 2 (d • V x d) + \ (fen - 2fci 3 ) (V • d) 2 + 

+ \k 2 2 (d-Vxdf + i (£; 33 + 2k 13 ) (d • Vd) 2 - 

- {h2 - k 2S ) (d • V x d) (V • d) + fci 3 V • (dV • d) - 

- \ (k 22 + M V • (dV • d - d • Vd) - 

- k 23 V • [(d • V x d) d] (11) 



where fe^'s are the elastic moduli (Frank's notation [5] is used in this paper), among 
which the relationship /C24 = \ ik\\ — ^22) holds. 

In relation (11), the first term originally represented the splay structure result- 
ing from the change of symmetry of the material, although liquid crystals having 
ferroelectric features are rarely found in nature. The elastic moduli k2, &12 and /C23 
are relevant to the chirality of the molecules. The three squared terms stand for, in 
sequence, the splay, twist and bending modes of the elastic deformation. The last 
three terms relate to boundary distortions. 

The viscous stress is supposed to be a linear function of the rate of strain tensor 
S and the rotational velocity of the director N [9,10] 

T v = ai (d • S • d) dd + a 2 dN + a 3 Nd + a 4 S + « 5 d (S • d) + a 6 (S • d) d (12) 

in which the independent scalar variables such as the mass density and the tem- 
perature have been absorbed in the viscosity coefficients. The viscosity coefficients 
satisfy the relationship [11] 

«2 + «3 = «5 — «6 (13) 
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The two considered variables are given by 



S = ^ (vV + Vv) (14) 

and 

N = d- x (v-r) x d (15) 

where r is the position vector, the overdot indicates the convected derivative. 
The movement of the director is governed by the equation [9,10] 



d_ 
~dt 




pddR = JJJ (G + g) dR + JJds-U (16) 



R R 



where p is the moment of inertia per unit volume (ML -1 ), G and g are respectively 
the external director body force and the intrinsic director body force, having the 
dimensions of torque per unit volume (ML _1 T~ 2 ), and II is the intrinsic stress 
(with the dimensions of ML -2 T -2 ). 

The magnetic field and electric field forces may, respectively, be obtained from 
the magnetization energy [9] 



W m = -^M H (17) 



and the polarization energy [29] 



W P = -^P-E (18) 



where P is the electric dipole moment. 

The constitutive relationship of P linear with the electric field is given as [9] 

P = ^— + — ~ e± (d ■ E) d = ^— + — (d ■ E) d (19) 
An An An An 

where ey and ej_ are the constant dielectric permeabilities with the interpretations 
indicated, and e a is the dielectric anisotropy. 

With the help of relations (6) and (17)-(19), the field-generated forces defined 
by [9] 

od 

are interpreted as 

G™ = - X a(d-H)H (21) 

G e = --!^(d-E)E (22) 
An 
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The intrinsic body force g consists also of three parts 

g = Ad + g e + g" (23) 

where A is a constant. 

The conservative and non-conservative body forces exerted on the director are 
[9] respectively 

dW d 

g ° = -5T (24) 

and 

g - = 7l N + 72 S-d (25) 

where 71 and 72 are the viscosity coefficients relevant to the intrinsic motion. 

Relationships between the 7 coefficients in Eq.(25) and the a coefficients in 
Eq.(12) are [9] 

71 = «3 - "2, 72 = «6 - «5 (26) 
The stress II has only a conservative part [9] 

,, ow d 

n = am (27) 

The extra-body force g e + g" and the extra-stress II generate a torque on the 
director. 

2.3. Energy balance 

The general form of the energy conservation law in the continuum is [9] 

!///fv.v + fa.a + ^* 

R 

'f ■ v + G ■ d + Q h ^j dR+ 

R ' 

+ jj cfe-^T-v + II-d - 3 h ^j (28) 



where U is the internal energy per unit mass, Q h is the heat supply per unit volume 
per unit time, and J h is the heat current flowing out of the volume per unit area 
per unit time. 

For uniaxial liquid crystals, the heat flux is supposed to have a linear relation- 
ship with the thermal field [9] 

J h = /3 ve + Pi (d • V0) d (29) 

where G is the temperature and (3$ and j3\ are the heat conductivities. 




6 



3. Adaptation to membranes 

Now we will adapt the Ericksen-Leslie dynamic theory to a bilayer membrane. 
Unless otherwise stated, the Latin suffix runs over 1, 2 and 3, whilst that the Greek 
suffixes run over 1 and 2. Italic symbols indicate the components of a vector or 
tensor in the moving local frame (see below) and block symbols represent those 
in fixed global coordinates (inertial reference system). Summation convention is 
assumed throughout. 

As suggested in the theory of thin shells, one defines a geometric surface situ- 
ated at the middle position in the membrane and names it the "middle surface" (see 
Fig 1). The word "surface" in the following means, uniquely, the middle surface. 

Let r be the vector from an origin in Euclidean space to a point in the middle 
surface. Any spatial position in the membrane is then determined by 

x 1 = r i (9 a ,t) + 9 3 n i (9 a ,t) (30) 

where x 1 denotes the spatial coordinates, 9 l are arbitrary parameters and n 1 is the 
normal vector to the middle surface. 

If 9 1 and 9 2 are identified with the curvilinear coordinates on the middle surface, 
then the local base vectors are given by 

dr e 1 x e 2 1 dr dr 

e 3 = n= - — n =^=t^tX— (31) 



d6 a |eixe 2 | T/E89 1 89 2 

where a is the determinant of the metric tensor 

a a p = e a -ep (32) 

The second- and the third-order magnitude of the surface are 

de 3 de a 
ba ( 3 = -^-e (3 = Mr e 3 (33) 



Ca(3 ~ 89* ' 89P (34) 



de 3 de 3 
89* ' 1)9? 

For uniquely determined surface, the fundamental magnitudes satisfy the iden- 
tity 

c Q/3 - 2Hb af 3 + Ka a p = (35) 

where 

2H = a^bap (36) 
is the first curvature (mean curvature) and 

K=p4 = - (37) 

\a a f3\ a 



the second curvature (gaussian curvature). 
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The director vector d (a unit vector parallel to the long molecular axes) is 
considered to radiate from the middle surface on the side to which the normal 
vector is positive (Fig. 1). 

The molecules displace in the Euclidean space at a velocity v. With respect to 
the local bases which move with the surface, the velocity has the components 

v a = v 3 = (38) 

because the particle flow is supposed to be parallel to the surface. With respect to 
the spatial coordinates x 1 the velocity is given by 

v 1 = i 1 = t>° + r 1 (39) 

where 

J - = w (40) 

is a time-dependent hybrid tensor. The tensor t l a is used to associate the flexible 
and movable surface with the inertial reference system. For example, if the director 
is known by its components d k with respect to the local bases, then its components 
in the global coordinates must be d k = t^d a + n k d s . Conversely, while one has the 
expansion of a field vector at a certain position of the middle surface with respect 
to the fixed coordinates, say, E k , one knows naturally, with respect to the moving 
frame, its tangential parts E a = t k E k and its normal part = E 3 = E k n k . The 
subscript (n) denotes the normal part of the vectors. 

One defines the varying configuration of the surface in a way that any point 
of the surface at the considered instant (^ 1 ,^ 2 ,t) is the endpoint of a translation 
of the point which occupied the same curvilinear coordinates on the surface at the 
preceding instant (6' 1 ,6' 2 ,t°) along the normal vector n°. The displacement speed 
of the surface is thus determined by 

where w is the rate of displacement. 

Across the thickness of the membrane the non-uniformities, such as that of the 
mass, the temperature and the internal energy, as well as that of the velocity and 
the stress, are supposed negligible. (The assumption v a 3 = implies that the two 
molecular layers do not mutually slip.) It follows that the differentiation of any 
tensor property of the material is a covariant differentiation along the surface. The 
surface derivative operator is defined by [30] 

V.=e*J, (42) 

Replacing the three-dimensional operator V in Eq.(ll) by the two-dimensional 
operator V s and multiplying both sides of the expression by the membrane thickness 
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h yields the surface density of elastic energy of the membrane. To the first-order 
terms (to small tilting angle of the director to the normal) , it is 



-k[ (^2H - d a d a - ^d a d a a^j + k' 2 e ap d f3 d ct + 
+ k' n H ^2H - d a d a - ^d a d a a^j - 



(k 22 + k 24 ) 



K-2H (d a d a + ^d a d a a^ + d aj p {T) b pc 



- 2k' 12 He al3 d (3 d ct - 2k' l3 d^d p H 



(43) 



where the symbol d a denotes the partial derivative with respect to the parameter 
6> a , e af3 is the two-parametric permutation tensor, the prime indicates the surface 
covariant derivative and the subscript (r) the tangential part (see Appendix A). 
Here the elastic moduli in Eq.(ll) have been replaced by 



k; — hkii 



k' 



hk 



t-3 



(44) 



Relationship (43) includes the energy contributed by the non-uniform orienta- 
tion of the molecules, which was emphasized to be important but excluded from 
consideration in the Helfrich's theory [7]. 

Letting d = n, one recovers Helfrich's formula [7] 



where 



Wf Q) = k c H + k' c H 2 + k c K 



k c — 2ki, k c — 2/c-q, k c — (k 22 + k 24 ) 



(45) 



(46) 



Hence the framework encompasses Helfrich's theory as a special case where the 
long molecular axes are parallel to the normal vector. Recently, Ou-Yang and Liu 
deduced Eq.(45) from formula (11) [31]. 

The first term in Eq.(43) and in Eq.(45) now stands for the non-symmetry of 
the two molecular leaves [7]. The last three terms in Eq.(43) and the last term in 
Eq.(45) are actually meaningful only when the membrane has a free edge such as 
a helical strips (e.g. [32,33]). Otherwise, the orientation of the molecules at the 
periphery is determined by the veins of the solid boundary on which the membrane 
is braced rather than the internal molecular field. For a closed vesicle, these terms 
are vacuous of physical content. 

The following kinematic theorem will frequently be used later, which is, in fact, 
an analogue of the Reynolds' transport theorem in two-dimensional fluids [28]. 

If <p is any function of position on a surface and of time, which can be any 
scalar or tensor component, and s is a material part of the surface, then 



ds 



(47) 
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where v a a + a/ (2a) is the dilation of area. 

The change in area of the membranes arises from, besides extension and com- 
pression, tilting of the molecules also [34-37] . For the surface configuration defined 
above, the rate of change of area is related to the displacement speed by [38] 

— = -2Hw (48) 

If there is neither a source nor a sink in the membrane, then 

v% = (49) 
On account of (48) and (49), one rewrites (47) as 



cf>ds= If (j)-2Hw(^ ds (50) 



d 
~dt 



3.1. Mass balance 

Let 7 be the surface density of mass. The mass conservation principle requires 

|//7* = (51) 

s 

Recalling theorem (50) and taking off the integration one obtains 

d t -f + i>%7 - 2Hw-f = (52) 

where the convected derivative d t = d t + v^dp has been taken into consideration. 

If is the mass percentage of the ith species inlaid in the membrane, then 
the balance of this species gives 

0t£(i) + v a d a £ (i) - 2Hwi {l) = Q {i) - Jg }> a (53) 

where is the surface chemical source and J(j) the surface diffusion flux of the 
species. 

3.2. Conservation of momentum 

The net force on an arbitrary area with periphery / resolved in the direction 
mp is given by 

j t J J \wds = J J F 3 ds + j> mpT 133 dl (54) 

s s I 

j t jj \v a ds = J J F a ds + j mpT^ dl (55) 
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Making use of theorem (50) and Eq.(52), one obtains Cauchy's equation for 
the displacement movement of the surface 

"fw = F 3 + T p3 (56) 

and that for the internal flow of the particles 

7 v a = F a + Tfj (57) 

In view of the perturbation analysis, one expands Eqs.(56) and (57) below on 
the assumption that the tilt angle of the director with respect to the normal vector 
is small. All the expressions in the previous section are valid once the spatial oper- 
ator V is replaced everywhere by the surface operator V s and the energy densities 
Eqs.(ll), (17) and (18) are multiplied by the membrane thickness. The details of 
the derivation are given in Appendix C. 

As is shown in Eq.(3), the body force F consists of the gravitational attraction 
and electromagnetic field-generated forces. The expansion of these vectors expressed 



in Eqs.(4)-(7) with respect to the local bases gives 

F g 3 = 19(n) (58) 

F g a = 7&4 (59) 

F m 3 = H Jjin J [ x± (H 1 - H^n 1 ) + XaH^tj^] (60) 

F m a = Hj.,14 [ X ± (H 1 - H^n 1 ) + x a H (n) t^d' 3 ] (61) 

F e 3 = le E {n) (62) 

F e a = 7e EitL (63) 
where 7 e is the surface density of electric charges. 



The stress tensor T consists of (as shown in (8)) an isotropic part, a conservative 
part and a non-conservative part. The application of the surface divergence operator 
to the isotropic stress yields 

T iP* = _ a ^ dpa (64) 

where a is identified with the surface tensor. The elastic restoring force of the 
membrane arises from the tilt of the molecules and the curvature of the sheet. 
When the long molecular axes have a small angle to the normal vector while the 
shape of the membrane is flat, then the surface force reads 

T e/ ^ = (65) 

T e % = -k' ia ^d a ^ {r) + k' 2 e a °a^d a , 7/J(T) (66) 

Conversely, when the membrane is flexible while the long molecular axes remain 
parallel to the normal vector, then the surface force is given by 

T e ^ = (k[ - 2Hk' n ) (4H 2 - 2K) + (k' 22 + k' 2A ) 2HK (67) 

T ep * = 2k[a al3 df3H - 2k' 2 e af3 d p H - 2k' u (b af3 + 2Ha af3 ) d p H+ 

+ (k' 22 + k' 2A ) a af3 d p K + 2k' 12 e aP {b a p d a H + d p H 2 ) (68) 
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One considers next the viscous stress. It is easy to verify that for the two- 
dimensional fluid the rate of strain given by Eq.(14) is 

S = -(v a ,(3 + Vf3, a -2wb af3 )e a eP (69) 
and the rotational velocity of the director given in Eq.(15) 

N : 



(70) 



where d a = d t d a + v a d a a ^ . In writing Eq.(69) one took into account the formula 
[38] 

d AM = -2wb\n (71) 

(see Appendix C.) Substituting expressions (69) and (70) into (12) gives the viscous 
resistance in a two-dimensional fluid. In the limiting case 2H — K = it is given 
by 

TV % (d„d° + ^d°d a a^j + v^ a{T) + v a d^ {T) - 



+ ^-E^ A/3 + V X , rf) d a + (V^ X + V X , M ) d a , p {T) ] 

T v%=p^(vx,^ + v,,xp) 
and in the other limiting case d || n it is given by 

T v ^ = -(jji + rj) AH 2 w + V [AKw + (b XtI - Ha x ») (v x , „ + a)] + 



(72) 
(73) 



rpv P& 

1 ,P 



+ f [v a ,f& + v% (2H)} 

- fia^dp (2Hw) + "Le^ (v x , + Vli , xp) - 

- 2 V [dpw {b ap - Ha ap ) + wa ap dpH] + 
+ a 2 Hv f3 b^ + ^ {2Hv p b pa - Kv a ) 



(74) 



(75) 



where \i is the dilatation viscosity, rj is the shear viscosity and E a ^ x ^ is the fourth- 
order tensor E af3x ^ = a aA a^ + a^a f3x - a ap a x ^. In deducing (72)-(75), formula 
(48) was used. 

The momentum balance for the intrinsic motion over the area s Eq.(16) is 
written as 



= JJ (G 3 + g 3 ) ds + j m^IT 33 dl 

s I 

^-jjz f d a = jj {G a +g a ) ds + j> m p U^ a dl 



(76) 
(77) 
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where 7 = 7/1. In the usual way, Eqs.(76) and (77) lead directly to the point 
expression of the law 



o = G 3 + < 7 3 + ir 93 



7(f 



f3a 



(78) 
(79) 



G a + g a + W J3 

When a polar fluid is exposed to the external magnetic field, the director is 
reoriented by the force G m , as given by Eq.(21) 



G m 3 = —XaH(n) (Hit^cT + #( n )) 

Cr n a = - X aH j t4(H 1 t^ + H (n) ) 



(80) 
(81) 



In an external electric field, it is subjected to the force G 6 , as given by Eq.(22) 



G e 3 = ~^ E (n) (Ei^cT + E {n) ) 
G e a = -^Ejt J Q {Eit l a (F + E (n )) 



(82) 
(83) 



For a two-dimensional fluid, the extra-body forces given by Eqs.(24) and (25), 
in the limiting case 2H = K = 0, are given by 



9 e3 = - k' 2 e^d p d a 



9 v3 =0 



g va =71 + \e xp d x v^d p e^ + ll E ^ (v\ jfl + v^ x ) dp 



and, in the other limiting case d || n, 

g e3 =o 



9 ea =k[ 3 a^d (2H) 
9 vS =0 
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(84) 
(85) 
(86) 
(87) 



(88) 
(89) 
(90) 
(91) 



where 73 is referred to the shear motion. 

The surface force to which the director is subjected during the elastic defor- 
mation, given by Eq.(27), takes the form 



=k' 2 e^d p d c 



(92) 



+ k' 22 e af3 dp (e^d x d^) + (k' 22 + k' 2A ) a^d" 



-k\ 
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e ap d p (d a d a + ^d a d a a}j + a af3 d (e XfI d x d^ 



(93) 



13 




for 2H = K = and 

= (k[ - 2Hk' n ) 2H + (k' 22 + k' 2A ) 2K (94) 
= - k'^a^dp (2H) + k' 12 e af3 dp (2H) (95) 

for d || n. 

3.3. Energy balance 

The energy balance Eq.(28) for any material area of the membrane reads 

Jt // 2 ^ Vct + w2 ) + \d a L+lUds 

F a v a + F 3 w + G a d a + Q h ^j ds+ 

+ j mp (T pa v a + T^w + U^da - J h ^ dl (96) 
i 

where Q h is the rate of local production of heat per unit area and J h is the surface 
diffusion current across a unit linear element of the periphery of the considered area. 

Taking account of Eqs.(56)-(57), and (78)-(79), one obtains from Eq.(96) a 
differential expression of the law 

1 U = Q h - J ha a + T^ a v a ,p + T^dpw + Tl^d^p - g a d a (97) 

For practical uses, one replaces the internal energy by a testable physical quan- 
tity and obtains 

7 c,e - T e/3c v - T eP3 dpw - n^da, p + g ea L 

= Q h -J ha a + ("^ a + T vPa ) v a , p + T v03 dpw - g va d a (98) 

where Ci is the heat capacity per unit mass in the static state [39], the convected 
derivative of the temperature is considered as 

Q = d t Q + v a tld j Q + wn i djQ (99) 

(see Appendix D.) 

The non- vanishing terms of Eq.(98) containing T e ' v , g e ' v and II are 

T ePa v a ,p = - (fcX ff d%(r) " k' 2 a^e^d^ x ) v a ,p (100) 
p = (k'y a + k' 2 e Pa ) d a , p (101) 
T vpa v a , p JlE*f»* (v x , n + v», a) v a , p (102) 



T v(33 dpw=\a 3 



d? + \ {e^d x v^ e^d a 
+ ^ E P*x» (uAj ^ + ^ a) da ^ d pw (103) 
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for 2H = K = 0, and 



(k[ - 2Hk' 11 ) b? a + (k' 2 - k' 12 ) b^ a + 



+ (k 22 + k'u)Ka (3a \v a ,p 
g ea d a =k' 13 d a a a % (2H) 
Ii pa d^ = \(k[ - 2Hk[ 1 ) a pa + (k' 2 - 2Hk' 12 ) £ /3q + 

+ (k' 22 + k' 24 ) (2Ha l3a -b f3a )]d a , p 
T v/3a v aj /3 =|(?7 — n) 2Hwa^+ 

+V \E a ^{v x ^ + v^ x )-2wb^ } Ua>/3 



T vf33 d p w = 



^v°b%w 



g va d a =^v p b^d a 
for d || n. Here the covariant derivative of d a is given by 

d a , f3 = d t (rf a ,/3(r)) + u ^a,7(r) + ^ 7 ^a, 7 /3(r) 

(see Appendix C.) 

For the local production of heat it is commonly written as [39] 



(104) 
(105) 

(106) 

(107) 

(108) 
(109) 



(110) 



(111) 



r k 



where Hk is the enthalpy of the kth species per unit mass, v^ r is the stoichiometric 
number of species k in chemical reaction r and b r is the reaction rate of chemical 
reaction r referred to a unit area. 

The divergence of the heat flux Eq.(29) is given for the two-dimensional fluid, 

by 

J h ° a = - «± (0,6) k (g jk - nm k ) + n a 2H (n^O) + 



+ K a 



t^O) (2Hd^ + bid-) - nJ ($0) k 4d°- 



- (n^G) (d a d a + ^d a d a a 



(112) 



where k± is the heat conductivity in the direction perpendicular to the director 
vector and K a = k\\ — k± the thermal anisotropy. 
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In the above derivations, one takes an arbitrary three-dimensional coordinate 
frame as the inertial reference system and two arbitrary lines on the middle surface 
as the curvilinear coordinates. For most purposes, however, one may choose a simple 
coordinate system. For instance, if one defines the two principal curvature lines as 
the curvilinear coordinates, which leads to an orthogonal conjugate system, then 
the expressions will be much simplified. 

4. Concluding remarks 

By adapting the electrohydrodynamic theory of uniaxial liquid crystals to thin 
films, one obtains a group of dynamic equations, including basic Eqs.(52)-(53), (56)- 
(57), (78)-(79) and (98), as well as relevant expressions (58)-(70), (72)-(75), (80)- 
(95), (100)-(109) and (111)-(112), which may conveniently be used in the study 
of hydrodynamic and electro dynamic phenomena of biomembrane matrix. The 
surface density of elastic energy Eq.(43) may be used to the hydrostatic studies of 
membranes with particular regard to molecule tilting. 
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Appendix A. Preliminaries 

The Gauss formula is given by 





and the derivative of the normal vector 



<9e 3 



(A2) 



dei 3 



where 




Some surface differential invariants of the director vector are given as 




V s d 
V s x d 



e a(3 [(<Vs + d*b a p) e a - df3d a e 3 ] 
e af3 (rf 3)/3 e Q - dpd a e 3 ) 




(A3) 
(A4) 



{Ah) 
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where e af3 is a two parametric permutation tensor, having the components e 11 = 
e 22 = 0, e 12 = — e 21 = 1/y/a. The subscript (r) distinguishes the tangential part 

d« p{T) = dpd« + d°{^ (A6) 

from the whole covariant derivative 

= dpd? + d° | ^ } - d% = d a p{T) - d% 

d s f3 =d^d 3 + d a b^ (A7) 

Supposing A is a field vector. The surface derivative of the vector is associated 
with the spatial derivative of the vector by 

V S A = A k ,/ a e Q (%eP + n k e 3 ) (A8) 

Let B be a second-order tensor (B = T or II). The contraction of the third- 
order tensor V S B gives 

V s B = (V; (r) - B&b% - 2M 3 j e a + 

+ (B^ (T) + B^b aP - 2e 33 ) e 3 (A9) 

with 

The normal vector to the middle surface has the spatial components 

Hi = \e^ a t^ {All) 

where is the three-parametric permutation tensor, which is equal to v /g when 
ijk is an even permutation of 1, 2, 3, equal to — y/g when it is an odd permutation 
and in other cases, g being the determinant of the metric tensor of the spatial 
coordinates gy. 

Appendix B. Surface density of elasticity energy 

Volume density of elasticity energy Eq.(ll) can otherwise be expressed as 

W d =h (V-d)-fc 2 (d-Vxd) + \k xl (V • d) 2 + 

+ \k 22 (d • V x d) 2 + 1 (fc 33 + 2k 13 ) (d • Vd) 2 - 

- \ (k 22 + M [(V • d) 2 - dV • dV 

- k 12 (V • d) (d • V x d) + 

+ A;i 3 d • V (V • d) - £; 23 d • V (d • V x d) (A12) 
17 



Replacing V by the surface operator V s defined in Eq.(42), and multiplying 
the moduli by the membrane thickness h, one obtains from Eq.(Al2) the surface 
density of the elasticity energy 

W d = - k[ (2H - cP a(r) ) + k' 2 e^d p d a + k' xl H (2H - cT q(t) ) - 

- (k' 22 + k' 24 ) (K - 2Hd a a[T) + d a ,ft {T) hP a ) - 

- k' 12 2He aP dpd a - k'wdPdp (2H) + O (<T 2 ) (A13) 
where k\ = hki, k'^ = hkij. 

Appendix C. Forces sustained in membranes 

CI. Forces being exerted on the director 

Cl.l. External forces 

For a two-dimensional fluid, the force to which the director is subjected to in 
an external magnetic field is given by 

G m = - Xa (d H)H 

= - Xa (H k £d* + H (n) ) (Hj^e + H (n) e 3 ) (AU) 

and that in an external electric field is given by 
G e = — ^(d-E)E 

= (E k tU a + E {n) ) (E/ Q e Q + E (n) e 3 ) (Alb) 
CI. 2. Extra body forces 

CI. 2.1 Elastic restoring force. Inserting Eq.(Al2) in relation (24), one obtains 
dW d 

= k 2 V x d - k 22 (d • V x d) (V x d) - 

- (fc 33 + 2k 13 ) d • (Vd) • (dV) + k 12 (V • d) (V x d) - 

- k 13 V (V • d) + k 23 [V (d • V x d) + d • V (V x d)] (A16) 

With respect to the local reference frame, it is expressed as 

g e =k' 2 e aP (d°b af3 e a - dpd a e 3 ) - 

- (k' 33 + 2k' l3 ) (2#6 Q % - Kd a ) e Q - 

- k' 12 2He af3 (d7b a pe a - d p d a e 3 ) + 

+ k' 13 a ap d (2H - d a d a - T^^a) e Q + 
+ k^ 3 a a % {e^dxdj e a + O (> 2 ) {All) 
18 



g e = 



Cl.2.2. Viscous resistance The rate of strain S given in Eq.(14) is written for a 
two-dimensional fluid 



S = - (vV s + V s v) 
= l - (fy W + e«a a v) 

= ^gi J (t J aV i )/ ,+tj 3 V i >Q ) 

= |gij k (*X + r 1 ) , + 4(^ + r i ). 



e Q e /3 



= ^ (fa, /3 + vp t a - 2w6 a/ 3) e Q e /3 



(A18) 



for [28] 



and [38] 



= -2^6 



A/j 



(71) 



The rotational velocity of the director N defined in Eq.(15) is written as 



N = d - - V s x (v - r) x d 
2 



e a + 



(A19) 



where 



d a =d t d a + v' T d% {Ty 
Inserting Eqs.(Al8) and (A19) into relation (25), one obtains 

g" =7i (d a - \v p b^ - l - £ ^d x v^dpe^ e Q - y 2 2Hwd a e a + 
(v x , n + \) - 2w (b^ - HaP") 
+ ^iv a b a pd p e^ + o(d a2 ) 



+ 72 



(A20) 



where, 7 2 is referred to the area dilation, 7 2 ' to the shear motion, the fourth-order 
isotropic tensor E a,3x ^ = a^a 13 ^ + a a ^a l3x — a XfJ, a al3 is used for the tensor operation 
compatible with the transverse isotropy of the material 
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CI. 3. Extra surface force 

The substitution of Eq.(Al2) into relation (27) gives 



n 



dW 



<9(dV) 

=k x I + k 2 d-e + k u (V • d) I - k 22 (d • V x d) d • e+ 

+ (k 33 + 2k 13 ) d • (Vd) d - 2 (k 22 + k 2A ) [(V • d) I - Vd] + 
+ fcia [(V • d) d • e - (d • V x d) I] - 

- k 23 [(V x d) d - d ■ (Vd) • e] (A21) 

where e is the third-order permutation tensor. With respect to the local bases it is 
given by 

II =& y 1 a /3a e /3 e Q + k' 2 e^ a (e^e a - d a epe 3 ) - 

- k' u ^2H - d a d a - ^dTdaoj a /3a e /3 e a - 

- k' 22 e x »d x d^ a epe a - (k' 33 + 2k' 13 ) d°1%epe 3 + 



+ 2 (k 22 + k 24 ) 



2Ha f3a - b pa - ( d a d a + Y a dadaa ' 



+a (3 ' y d t 



7(r) 



e^e a + & /3a d Q e /3 e 3 j- 



12 



2ife /3c 



+^<9 A d M a^ e /3 e Q + 2#£ Q ^ Q e /3 e 3 }- 
- /^e^teeg + O (ef* 2 ) 
The application of the surface divergence operator to Eq.(A22) yields 



(A22) 



=k' 2 e /3rT d a b% - k'^a^dp ( 2H - d n d a - —d a d n a ) - 



1 

2a 



- k' 22 e^d p (e^dxdj + (k' 33 + 2k' 13 ) (2Hd% - Kd a ) + 

+ 2 (k' 22 + k' 2A ) Kd a - 2Hb f3a dp - aP a d p (d^ + ^d a d a a^ + 

+a^d c 



,7/3(r) 



k' 



12 



1 



e^dp {2H - cU CT - 7^d°d a aj + 
+a /3a <9 /3 (e Xfl d x d^ + 2He pCT b a p d (T + O (V 2 ) 



(A23) 
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Il^g =k[2H - k' 2 e pa d p d a - k' n 



AH 2 - 2H ( d a d a + Y a d ° d(ja 



- (k' 33 + 2k' 13 ) \bZ<r tfKr) +d%(2H) 



+ 4 (k' 22 + k' 2A ) 



+ 



K-H( d a d a + ^d a d a a ) + 



+df 3 dpH + b< 3 °d a ^ {T) 
+ k' 12 e^d a df3 (2H) + k' 23 2He pa d p d a + O (V 2 ) 

C2. Forces exerting on the continuum 
C2.1. Gravitational attraction 

The gravity given in Eq.(4) is written for the two-dimensional fluid 

F 9 = 7 g = 7gk {t k a e a + n k e 3 ) 



(A24) 



(A25) 



C2.2. Ponderomotive forces 

The force the external magnetic field exerts on the body of the polar fluid, 
given by Eq.(5), takes the form 

F m = H kJ [ x± (ff - njH (n) ) + XaH (n) 4^] (£e a + n k e 3 ) (A26) 
In deducing (A26), one used relation (6) and the formula [28] 

tj^ = g« - nV (A27) 
The Lorentz force in an electric field (7) is given by 

F e = 7e E = 7e (E k t k a e a + E {n) e 3 ) (A28) 



C2.3. Isotropic force 

The application of the surface divergence operator to the isotropic stress ex- 
pressed in Eq.(9) gives 



(A29) 



C2.4. Elastic restoring force 

Inserting Eq.(Al2) into relation (10) yields 

dW d 

T e = - (Vd) 



<9(dV) 

- fciVd + k 2 (Vd) x d - k n (V • d) (Vd) - 

- k 22 (d • V x d) (Vd) x d- 

- (A; 33 + 2k 13 ) (Vd) • (dV) • (dd) + 

+ 2 (k 22 + k 2A ) [(V • d) Vd - (Vd) • (Vd)] + 
+ k 12 [(V • d) (Vd) x d + (d • V x d) Vd] + 
+ k 23 [(Vd) x (dV) • d + (Vd) • (V x d) d] 



(A30) 



21 



With respect to the local bases it is given by 
T e =k[ [(V Q - a^cP 7(T) ) e/3 e Q - d a b^epe 3 



+ 



+ k 2 e 

- k' 



/ _cra 



2#6 /3a - -^=d a (d a y/a) b pa - 2HaP a d a a{r) 

-2Hb^d a epe 3 \ + k' 22 e^d x d^ a<J b^epe a - 
- (k' 33 + 2k' 13 ) (2Hb f3a d a - Kd 13 ) e p e 3 + 
+ 2 (k' 22 + k' 2A ) | Ka^ + a^d aMr) b aa - -±=d a (<Ty/a) 

+ (b^ - 2Ha^) d% (T) epe a - Kd%e 3 |- 



+e x »d x d^ 



2He^b P x d ll epe 3 \- 



- k' 23 Ke^d a epe 3 + O (> 2 ) (A31) 
The dot product of the surface derivative operator and the stress T e gives 
T e/ ^ =k[ ^ a dp (2H) - d a af3{T) a pa + 2Hb f3a dp - Kef*] - 

- fc£ {e^ [dp (2H) + Kdp - a^d Pna(T) ] - 2He^by^} - 

- k' u | (bP a + 2Ha f3a ) dp (2H) - a^d^^dp (2H) - 

- (b pa dp + 2a pa dpH) 

+2H (2Hb pa dp - aP a d% m - K<P} | + 

+ k' 22 [s a °b% (e x *d x d») + e^dp (2H) {e^d x d^)] + 
+ (k' 33 + 2k[ 3 ) {AH 2 b pa dp - 2HKd a - Kb ap dp) + 

+ 2 (k' 22 + k' 24 ) {a^dpK + Kb^dp + {b^d^^) ^ + 
+ (b^-2Ha^) d« a p {T) - 



v ° 



+ 



(b ap dp + 2a ap dpH) 



d a (^d a ) 



- k' 12 \2e f3a (b a p d a H + dpH 2 ) + AH 2 e p(7 byp- 
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,/3a 



(2dpH + b%d a ) (±=d 1 (^ 7 )) + 



/3,<J7(r) 



+2a ia d 1 Hd f3ja(T) + 2Ha" T d 
- (b^dp + 2a^d p H) [e^dxd^) \ + 



2HKdf- 



(A32) 



T e ^=k[ AH 2 -2K -2bPd° f)(T) -d a d a (2H) + k' 2 e pa dpd a (2H) - 



-k' 



ii 



{AH 2 - 2K) ^2H - d a d a - ^-d° r <9 CT aj - 
- (d*&g + 2HdP) dp (2H) - AHb^p {T) 



(k 33 + 2k 13 ) 



2Hb^ p{T) + (&gd* + 2Hd?) dp (2H) 



+ 2 (k 22 + k 24 ) 



-d p dpK - K (d a d a + 7^% a ) 
2HK - d p dpK + 2Hb^d ajP(T) - 
- (AH 2 + K) (d a d a + ^d a d a a^j 
+ 2K[Je^ [(b°d a H + d a H 2 ) dp + Hbid^p {T) ] 



+ 



-(2H 2 -K) e^d x d^ + 
+ k' 23 e af3 (d a dpK + Kd a ,p {T) ) + O (> 2 ) 



(ASS) 



C2.5. Viscous resistance 

With the help of Eqs.(Al8)-(Al9), the viscous stress Eq.(12) is written as 

T v = - ^Hwa af3 + 2r)w (b ap - Ha ap ) - ^E a/3X ^ (v x , „ + V \) + 

+ (a 3 cP - l -v% - l - e^dp 

l -E a ^ (v x , „ + V x) dp - 2w (b a % - Hd a ) J e Q e 3 + 



- a' R 2Hwd a + 



23 



+ 1 «2 



al2Hwd l3 + 



+ 



X -E^ (v x , „ + V A ) rf a - 2™ (d a &«0 - d*) 



| e 3 e^+ 



(A34) 



where the viscosity coefficient «4 in Eq.(12) has been split into the dilatation vis- 
cosity \i and the shear viscosity 77, the a' coefficients refer to the dilatation of area 
and the a" coefficients to the shear motion. 
The dot product of V s and T" gives 



V*" = - fia^dp (2Hw) + V -E^ (v x , ^ + VfJb , Xf3 ) - 
- 2 V [dpw (b pa - Ha pa ) + wa 
+ a 2 Hv a b aa + ^ (2Hb aa v a - Kv a ) - 



«2 



- QJ 3 



d% + \bP a (v a d?) g (r) + d a v°d a H- 



+ 



+ a' 5 4H 2 wd a + a' 6 2Hwb al3 dp- 

- a'i [HE"*"* (v x , „ + x) dp - 4Hw (dpb? a - Hd a )] - 
l E a ^ (v x , „ + V a) d a b% - 2w (Hb a °d a - Kd a ) 



— a 



+ 



+ d 



(A35) 



iv /33 



= - iaAH 2 w + r] (b x » - Ha x ») (v x , M + v„ t x) - 
- 4 V w (H 2 -K)- — (b%v a jl3 + 2v p dpH) + 

+ a^3 t (dpd? + ^d%a)j + v°pd% {T) + 

+v' T d P af3{T) - X -e^ (d a dp + dpd a ) (e x »dxv„) }- 



d p d p (2Hw) + 2Hw (dpd 13 + ^-d^dpa] 



2a 



+ 



+ a '^E^ [(v,,x + vx,,)d a ]^ {T) - 
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a,f3(r) 



-20/3W (d a b af3 - Hd 13 ) - 2w 

+d (3 dpH - H (dpd + ^fdpa^j j 
(a 2 + « 3 ) Hv a b aP d p + O (> 2 ) 



b /3a + 



(A36) 



Appendix D. Energy balance 



In the following, one looks for a testable parameter to alter the internal energy 
U in the energy balance equation. 

With the help of Eqs.(8), (9) and (23), Eq.(97) is written as 

7 f/ _ (T e/3c X p + T'&dpw + U^d a , p - g ea d a ) 

= Q h -J h " a + [(-aa^ + T^ Q ) v a , p + T v03 dpw - g va d a ] (A37) 

From the first law of thermodynamics, one knows that 

dU-J2 ^ d h = dQ + dW dlss (A38) 

i 

where Wdiss is the dissipated work, k the ith work coordinate, Li the conjugate 
work coefficient, and U and Q are defined as usual. 

As U is the characteristic function of the temperature 6 and of the work 
coefficient k : one has 

Hl)-£(f) ei *. 

Equality (A38) may thus be written as 



^1 



dh. 



dQ + dW dl 



(A39) 



From the fundamental equations 

ci£/ = GdS + L idh 

i 

dF = -SdO + J2 L idh, 

i 

one has 



dS 



"1/ i 



-) =e, 

dh J t . \dl 
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where F is the Helmholtz function and S is the entropy. Substitution of the two 
expressions above into Eq.(A39) leads to 

dde-eJ2{^) dh = dQ + dW diss (AAO) 

where 

is defined as the heat capacity at constant work coordinates [39]. Dividing Eq.(A40) 
by the infinitesimal time interval, one obtains the instantaneous energy balance 



dt ^ \dO J, dt dt dt 

We turn back now to Eq.(A37). The terms in the parentheses on the left-side of 
Eq.(A37), i.e. the rate of change of the work of the elastic deformation, is identical 
to the second term on the left-side of Eq.(A42), whereas the right-sides of Eqs.(A37) 
and (A42) are identical. It follows thus that 

iQQ - T e/3 % Q , p - T^dpw - U^da, p + g ea d a 

= Q h -J ha a + (~aa^ + T^ a ) v a , p + T^dpw - g va d a (A43) 
Here the convected derivative of the temperature is given by 

9 = d t Q + (v a tl + «m j ) djO. 

The internal energy U in Eq.(A37) has now been replaced by the product C;6. 
The heat flux Eq.(29) is expanded with respect to the local bases 

J h = - djQ { (n ± tl + K a iiU a ) e a + ( Ka tld a + K||n j ) e 3 } + O (> 2 ) (AAA) 

where the coefficient (3q has been replaced by the heat conductivity perpendicular 
to the director vector — k± and the coefficient j3\ by the thermal anisotropy — K a = 

K±_ - K\\. 

The divergence of the heat flux is given by 

J ha a =J h : a(T) -2Hj» 3 

= - «L (6>je) jk (g kj - n k n j ) + K a 2H (n^jO) + 



- (n^jG) (d a + Y a d «°) 



-n J (5j6), t k d Q 



+ O (d° 2 ) (AA5) 
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Figure captions 

Fig.l Illustration of the middle surface and the local coordinate axes. 
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